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Abstract 

We compute a novel type of large non-Gaussianity due to small periodic features in general 
single field inflationary models. We show that the non-Bunch-Davies vacuum component 
generated by features, although has a very small amplitude, can have significant impact 
on the non-Gaussianity. Three mechanisms are turned on simultaneously in such models, 
namely the resonant effect, non-Bunch-Davies vacuum and higher derivative kinetic terms, 
resulting in a bispectrum with distinctive shapes and running. The size can be equal to 
or larger than that previously found in each single mechanism. Our full results, including 
the resonant and folded resonant non-Gaussianities, give the leading order bispectra due to 
general periodic features in general single field inflation. 



1 Introduction 



Recent progress in the theory of inflationary density perturbations has revealed a variety of 
mechanisms that can be responsible for large non-Gaussianities [1]. They are discovered not 
only in novel models, but also in simple models where such mechanisms were not previously 
anticipated. For example, while it is certainly expected that small periodic features in an 
inflation model can give rise to small oscillatory corrections to the power spectrum, it is 
recently found that the same features can boost the amplitude of the non-Gaussianity by 
several orders of magnitude [2]-|6] . These features introduce small oscillatory components to 
the inflationary background. The inflaton quantum modes can resonate with this background 
while they are still within the horizon, generating large interactions. These models are called 
the resonance models, and the resulting non-Gaussianity has a very distinct form. In this 
paper we give another example in which novel large non-Gaussianity can be generated due 
to periodic features. 

Features generically perturb the inflaton away from its Bunch-Davies vacuum state. The 
mode function of the inflaton is generally a linear superposition of a positive and a negative 
energy component. The usual Bunch-Davies vacuum is chosen to be the positive one because 
it asymptotes to the vacuum state of the Minkowski spacetime in the small wave-length 
limit. However any disturbance due to features, sharp or periodic, generically introduce the 
other component [SlH]. This component has to be small to be consistent with the observed 
approximately scale-invariant power spectrum. Therefore it can be neglected in the previous 
computation of the resonant non-Gaussianity. 

On the other hand, the negative energy component can enhance the folded bispectra [7], 
or more generally the higher correlation functions containing folded sub-triangles in the mo- 
mentum configuration [H]lll With such a component, the phase of the mode function product 
becomes stationary for the folded momentum triangle, producing a large non-Gaussianity. 
In Ref. [THTT]. the simple forms of non-Bunch-Davies vacua are put in by hand, largely 
motivated by the discussions of the hypothetical trans-Planckian physics [12]. For example, 
the inflaton wave-function at the transition point from Bunch-Davies to non-Bunch-Davies 
vacuum does not solve the equation of motion. Feature models provide solid examples of 
such vacua with self-consistent forms. In the resonance models we consider in this paper, 
although the amplitude of the non-Bunch-Davies vacuum is very small, it is generated deep 
within the horizon. As we will see, this dramatically enhances the folded non-Gaussianity 
through the interaction terms in the action that are least suppressed by slow-variation pa- 
rameters. These interaction terms can be different from those important for the resonant 
non- G aussianity. 

We comment that the type of non-Gaussianity we compute here arises not because we 
introduce new ingredients in the model. It exists in known types of models but previously 

"'^Thc folded bispectrum means the bispectrum evaluated at the folded triangle configurations (fci + fc2 — 
/ca = or cyclic). In higher correlation functions, the folded sub-triangle configurations are defined to be 
the configurations in which three momenta ki, kj and kij = |k,: + k^l form a folded triangle, for any i, j. 
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unnoticed. This was pointed out in Ref. [T]. The effect will turn out to be most dramatic 
in the framework of general single field inflation models [7]. So we consider the inflaton 
Lagrangian as a general function of the inflaton </> and its first derivative [T3] . 



£ = P(X,0), (1.1) 

where X = —^g^'''d^(j)du(t>- Small periodic features can come from different places, for exam- 
ple, the potential for slow-roll inflation [3lll] or the warp factor in the kinetic term for DBl 
inflation 0. To be general, irrespective of the mechanisms, the only quantities that we will 
keep track of are its oscillatory frequency u and some general oscillation amplitudes that we 
will parameterize later. Other notations are the same as in [7|. The background metric is 

ds'^ = -dt^ + a^d^^ . (1.2) 

The three slow-variation parameters are 

The two parameters that describe the higher derivative kinetic terms are 

2 _ P,x A _ 3XPxx + 2X^Pxxx ,n 

~ Px + 2XPxx ' S ~ 3Px + QX^Pxx ' ^ ' ' 

where subscript "X" denotes the derivative with respect to X. Both and A/S are kept 
arbitrary. The scalar perturbation ( comes in when we perturb the scale factor as a^e^'' in 
the comoving gauge in which the scalar field has no fluctuation. It is conserved after horizon 
exit and is related to the CMB temperature fluctuations by C ~ — 5AT/T. We will calculate 
its correlation functions. We will write the bispectrum in terms of a function S{ki, k2, k^) 
defined as follows [1], 

(C^) ^ Sik^,k,,ks) J P^i2ny6\J2^.) , (1.5) 

^ ' i=l 

where = //^/(Svr^eCs) is the leading order power spectrum. 
2 Resonant non-Gaussianity 

We first compute the resonant non-Gaussianity ^2j in the present context of the general single 
field models. For this purpose, the most important terms in the cubic action of the scalar 
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perturbation [7] are 



5*9 D / dtd^x 



a^e d 



ci 



2A 



S J 



(2.1) 



This is because these terms contain either the highest time derivatives of the background 
parameters H and Cg, or the least number of factors of a (in terms of conformal time dr = dt/a 
and C' = o-C)? both contribute factors of u/H. The first and third term are more important 
if the oscillation originates from e, while the second and third term are more important if it 
is from Cg or A/S. This gives the following three-point function, 




dr 



a^e d f T] 
c? dt \ c? 



Vk2^ + — ^(k2 ■ ka) ||Mfc^(r) 



2A\ -i-r dul 

n 



dr 



(27r)^5(^ ki) + 2 perm. + c.c. 



(2.2) 



where the "2 perm." denotes two other terms with permutation of indices 1, 2 and 3. The 
first term is a generalization of that used in [2] . But in general single field infiation, the other 
two terms can also be the leading terms. The dominant component in the mode function Uk 
is the positive energy mode 



iH 



Uk\j) 



[1 + ikcsT)e' 



-ikcsT 



(2.3) 



We separate the coupling into a non-oscillatory part and an oscillatory part, and denote the 
amplitude of the latter with a subscript A. For example. 



:s(t)"l(t)/l(|)/'"''^" 



(2.4) 



and similarly for the others. We assume that all the oscillatory parameters have the same 
frequency and phase, otherwise simple modifications are necessary, especially for the phase. 

The case with non-oscillatory couplings and Bunch-Davies vacuum (12. 3p in general single 
field models is considered in [7]. In this section we focus on the oscillatory part, and in the 
next section we study the effect of the non-Bunch-Davies correction. 

In resonance models, a common integral that we encounter when we compute the non- 
Gaussianity or wave function of mode function is of the form 



dre 



fir) 



(2.5) 



This integral can be interpreted as follows ^ . The first oscillatory factor is from the couplings 
such as (12.41) and has the background frequency. The second oscillatory factor is from 
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the mode functions when they are inside the horizon where is a comoving momentum. 
Due to the background periodicity, the phase of the integral repeats after 2nH/uj e-fold 
{ANe = -HAt, t -i-t + 2ti/uj, t -> re'^''^/'^ and K Ke^'^^'^). So the integral must be 
proportional to exp(±i(c(;/iJ) Ini^), where the plus/minus sign comes from the fact that a 
larger K resonates later with a larger t, corresponding to a larger /smaller phase as we can 
see from (12. 5p . The physical frequency of a i^-mode decreases as it is stretched by inflation, 
and the resonance happens when it coincides with u. So the resonance point r^, is given 
by d{KcsT) / dt = u, i.e. t^, = —uj/{KcsH). The leading contribution of the integration is 
obtained when each K-mode sweeps through this resonance point. We approximate it as 
a step-like function at with a width and an amplitude that we estimate in the following 
way. Once the mode frequency differs from u by Au, the integral cancels out if performed 
over Ati ^ 27t/u. On the other hand, it takes At2 ~ Au/{uH) to stretch a mode until 
its frequency changes by Au. Equating Ati and At2 gives the duration At ^ —Hr^^At ~ 
\fZjjH /[Kcg)., namely the resonance width. It also gives the number of resonance periods, 
u At/ {271) ^ y^uj/{2'KH). Since one period of resonance gives 2tt/{Kcs), the final resonance 
amplitude is ^j2/mjjjHl (Kcg). The rest of the integrand, denoted as /(t), varies much slower 
within the resonance width and its contribution can be simply approximated as /(t*). To 
summarize, the integral (12. 5p can be approximated as, up to a i^-independent phase ipo, 



fjr*) 2-KUj 
Kc. V H 



Elf \ Kc J —{t + 

\l 2u KcM 




(2.6) 



where Erf is the error function and used to model the step-like behavior. The above in- 
terpretation becomes the stationary phase approximation when we precisely integrate (12. 5p 
using t = —H~^ln{—HT) [5]. This gives the precise numerical numbers above and the 
phase ipo = ±{u/H){l — ln{uj/cs)) ± 7r/4. It also shows a more complicated behavior with 
superimposed small oscillations. We will use (12. 6p as the leading order approximation. 

Large resonance happens only if u ^ H. The terms with the highest order of /cjC^r 
dominate in (12. 2p . Other terms can be larger in the very squeezed limit ki/kj < H/u, but 
they are sub leading terms because their contribution to /atl's are suppressed by 0{H/u). 
We neglect these terms here. Using (12.60 in (12. 2p . we get the following bispectrum. 



5res = /]vtsin(|:lnir + (^,es)+/. 



ros2 
NL 



cos 



i-res3 ^1^2^3 . f ^ ^ , 
+ Inl sm [^hlK + V, 



(2.7) 
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where K = ki + k2 + k^, (fir 



-7r/4 + {co/H){l - ln(w/c,)) and 



^resl 
J NL 



fves2 
J NL 



a/tt 1 d f 7] 



/ W \ 1/2 



8^2 Hdt \cy A 
Vvr sa 



res3 
NL 



3V? 



4^2 



- 1 



2A 



a;\5/2 



(2.8) 
(2.9) 
(2.10) 



The first component is a generalization of |2] and has the same running behavior. The second 
has an additional mild shape. The third has an additional equilateral shape. All of them 
can be important depending on the parameters. 



3 Folded resonant non-Gaussianity 

Features also generate a non-Bunch-Davies vacuum component, which creates small oscilla- 
tory features in the power spectrum. In the resonance model, long after the mode resonates 
with the background, the background frequency is too large to have any further significant 
impact on this mode. Therefore, as in the sharp feature case, we can write the mode function 
as a linear superposition of two modes [5], 

?ifc(r) = u+{t) + c_M_(r) , (3.11) 

where m+ is the leading behavior fl2.3p and is the negative energy component, 

u_{r) = -j^={l - tkcsT)e''''^^ . (3.12) 

The coefficient c_ is small and its time-dependence can be solved using the linear equation 
of motion, 

v'l + cl^Vk - = , (3.13) 
where Vk = zuk, z = a\/2e/cs, and 

2a^H^l + 6) . (3.14) 

z 

The S denotes 

e 3r7 3s en es rf tis f] s . 

6 = h — ^- h— ^ , (3.15) 

24 2 428 2 2 AH 2H ' ^ ^ 



z" 
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in which the last two terms are usually the most important for our purpose. We denote 



6 = 60 + 6a sin(u;t) 



(3.16) 



where the second term is the fast oscillatory component and Sa denotes its amplitude. We 
assume that the oscillation frequency in fl3.16p is the same as that in (12.41) . as usually 
happens. We also ignored possible constant phases in (I3.16p . Otherwise it is simple to 
modify correspondingly. The integration we get from 



dr 



26aH 



UOT 



-2ikc,T 



cos{cot) 



(3.17) 



is of the same type as in (12. 5p . and the m_ component remains approximately zero until it 
gets excited at the resonance point [HIS]- Using the approximation (12. 6p . we get 



2tcoa — exp 



Ink + ip. 



Erf kc. 



2H , u 

(rH , 

00 ^ 2kc,H 



+ 1 



(3.18) 



The amplitude of the negative energy mode should be small. So it is negligible in the 
computation of the resonant bispectrum. However a larger resonance also makes the mode 
generated deeper within the horizon, at a scale 2H/uj smaller than the horizon. The deeper 
into the horizon the negative energy component is created the larger the associated folded 
bispectra [7]. 

Let us look at the details for the case where 1/cj. — 1 or A/S is equal to or greater than 
0{1). For the purpose of this section, the three most important terms in the cubic action in 
Ref. \7] are 



D / dtd^x 



lid 



1 



2A 



(3.19) 



because they are least suppressed by the slow-variation parameters. 

We use the first term as an example. We first consider the effect of the cutoff introduced 
by the step-like function in (I3.18p . Using m_ in one of the mode functions and considering 
the leading non-oscillatory part of the coupling, the oscillation factor in the mode product 
is cancelled for the folded momentum triangle. Since the U- component is generated at 
a scale ~ H/uj times smaller than the horizon size and the integral is proportional to r^, 
the enhancement is ~ [u/H)^. Multiplying the small amplitude in c_ and considering the 
pre-factor (l/c^ - 1 - 2A/S), we get f^i ~ (1/c^ - 1 - 2X/E)6a{uj/H)'^/\ 

The finite width in the step-like function in (I3.18P is also important. It makes the shape 
of the bispectrum sharply concentrated near the folded limit. To see this, we note that, 
away from the folded limit, the mode function product oscillates with a period Arphasc ~ 



-ki + k2 + ks 



or its cyclic. If this period is much larger than the width of the step- 
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like function, Ar ~ uo / {2H) / {kiCs) ., the component can be approximated as being 
turned on by an infinitely sharp step-function, and the integration around the sharp edge 
contributes significantly to the bispectrum. But in the opposite limit, the u_ component 
is turned on very slowly, and the fast oscillation from the mode functions averages out the 
effect of the step. Therefore, if we use the quantity —1 + k2lki + k^/ki (and its cyclic) 
to describe the deviation from the folded triangle limit, the bispectrum is suppressed if 
— 1 + A;2/^i + ^3/^1 > \/H/liJ (and its cyclic). 

The following is the full leading order results for fl3.19p . The bispectrum 5* is a summation 
of two terms, S\ and S'c, 



Sr. 



3V2n ( 1 



32 

k2kz 



2A 



kl 



16 



a[x\) cos(— In k\ + (p_) + h[x\) sin(— In k\ + 
n H 



+ 2 perm. , (3.20) 



coef(«,j,fc)-p- 



HJ 

c{xi) cos(— In ki + (/?_) + d{xi) sin(-^ In ki + (f_) ,(3.21) 
H H 



where coef{i,j,k) = 3/2 {j ^ k ^ i), —3/2 (j 7^ k,j{oT k) = i), 1/2 (j = k), and we sum 
over all 27 permutations. We have also defined xi = {1 — k2/ki — kz/ki){uj /2H) and cyclic, 
and 



a[x] 



c[x] 



1 



X'' 



2 + ix' 



H o, 
2 x^) 



e 2" COS X 



[X 



H 



H^2 



2 - —x^-^ --^^ 



(2x H x-^) 

CO 

2H 



e 21^ sm X 



e 21^-" sinx + (2xH x )e 

CO 



x"^ 
1 



e 21^ cos X 



xe 21^ sm X 



.-H_ 2 . _H_ 2 

-e smx + xe 2^^ cosx 



_H_^2 



cos a; 



(3.22) 
(3.23) 

(3.24) 
(3.25) 



These functions are the leading terms for 00 /H ^ 1. Note that fl3.22p - fl3.25p are all finite at 
the folded limit x — )■ 0. The phase ip- is related to ^^res in (12. 70 by ip- = ipres + i^/H) ln2 -|- 
7r/2. For Sc, we integrated over the terms with the highest order of kiCsT. The other terms 
are the subleading contribution which can become important in the very squeezed limit 
{ki/kj < H/uo). It is easy to see that this bispectrum is periodic-scale-invariant, under 
exp{2TTnH/co) where n is an integer, due to the symmetry of the Lagrangian. 
Since both the shape and running of this bispectrum are important, instead of plotting 
S{1, k2, k^) as usual, we plot the shape and running separately. The shape is the dependence 
on the momenta ratio X2 = k2/ki and X3 = k^/ki with fixed K = ki + k2 + k^. So we plot 
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1.0 °° 



Figure 1: Shape and running of a folded resonant bispectrum Sx'- the 2nd hne of f l3.20p with 
oj/H = 100, if - = and K = 3 (for the shape). 




Figure 2: Shape and running of a folded resonant bispectrum S^- the 2nd line of fl3.2ip with 
uj/H = 100, = TT and K = 3 (for the shape). 

S{K/ {1 + X2 + X3) , KX2I (1 + X2 H-Xs), Kx^/ {1 + X2 + x^)). The running is the dependence on 
K with a fixed shape which we choose to be ki = 2k2 = 2k^. So we plot S{K/2, K/i, K/4). 
The folded shape and the resonant running that we have discussed qualitatively can be 
clearly seen in Fig. [1] & |2l and are nontrivially combined. This will become clearer after we 
describe it with a simple ansatz. 

4 Simple ansatz 

With the upcoming high precision CMB data, either to test the Gaussianity of the primordial 
density perturbations or to search for possible large non-Gaussianities, it will be effective to 
perform general search schemes [TiHTT] followed by targeted optimal analyses P^VET] . Simple 
analytical ansatz of bispectra are necessary for these purposes. The bispectrum we find here 
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1.0 °° 



Figure 3: Shape and running of the ansatz (14.261) with uo / H = 100, = and K = (for 
the shape). 



is orthogonal to the scale- invariant folded bispectrum ansatz, and it is also approximately 
orthogonal to the resonant bispectrum ansatz (the first term in (12. 7p ) due to the leading order 
non-trivial shape. So a new ansatz is necessary. We find the most important properties of 
the bispectrum (13.201) and (13.211) are captured in the following simple form. 



-•fold-res 

Jnl ^-^P 
2 perm. , 



^3/5 



4 



ki ki 



sm 



C 
2" 



-1 + 



ki ki 



+ 2\nki] + (f 



(4.26) 



where C = u/H and the size of f^^^~^'^^ is given by the first lines of (I3.20p and (I3.2ip 
with proper normalization. An example is shown in Fig. [31 The exponential factor has 
been chosen to fit the particular decay behavior in ( I3.22p - (]3.25p away from the folded limit, 
namely, ~l/xasl<a;< \fujB. and ~ 1/a;^ (or 1/a;^) as a; > \Ju/ H, and may be 
adjusted. Such a bispectrum may be constrained or detected using the method of mode 
decomposition [TH4T6] by properly choosing the set of base modes. 



5 Size of bispectrum 

From (I3.18p . we know that the correction to the power spectrum is 
5PJP^ = 2V2^6A{H/ujf/\os{{uj/H) \nk + ^_) 



(5.27) 



The periodic oscillation can originate from that in H or Cg, and the time derivatives on 
these parameters contribute factors of u/H. So to keep 5Pc_/Pc^ small, we need to keep 
{Ae/e){uj / HY^'^ and {Acs/cs){uj / HY^'^ small, where Ae and Acg denote the oscillation am- 
plitudes. 
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On the other hand, the amphtudes of the resonant bispectrum and the folded resonant 
bispectrum are given by 



Here 1/c^ stands for either 1/c^ — 1 or A/S, and Acs/cs stands for their relative oscillation 
amplitudes. Note that l/c^ — 1 and A/S do not have to be much larger than one. They 
can be of order one or somewhat smaller than one, but the /atl's are still large. As we can 
see, depending on parameters, the folded resonant bispectrum can be the dominant one in 
the theory, equally important as or more important than the resonant bispectrum (12. 7p and 
the equilateral bispectrum in [7] . For slow- roll inflation, the terms that we computed for the 
folded resonant non-Gaussianity are no longer the (only) leading terms. For example, for 
models with only the canonical kinetic term, a similar calculation with the other terms in 
the cubic action shows that fj^j^'^'^^ ~ e5A{(^/ HY^"^ ~ Ae{co / H)^^"^ . Comparing to the scale- 
invariant part J'nl ~ the folded resonant non-Gaussianity can still have an enhanced 
amplitude, although in this case the resonant non-Gaussianity will be the most dominant 
component, fj^l ~ {At/ e){uj / Hf/^ . 

In summary, the overall results of this paper, the resonant non-Gaussianity (12. 7p and the 
folded resonant non-Gaussianity (I3.20p and (I3.2ip . give the leading order bispectra due to 
general periodic features in general single field inflation. 
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